Was the Universe Reionized at Redshift 10? 



Abraham Loeb 

Astronomy Department, Harvard University, 60 Garden Street, Cambridge, MA 02138; 

aloeb@cf a. harvard, edu 

Rennan Barkana 

School of Physics and Astronomy, The Raymond and Beverly Sackler Faculty of Exact Sciences, 
Tel Aviv University, Tel Aviv 69978, ISRAEL; barkana@wise.tau.ac.il 

Lars Hernquist 

Astronomy Department, Harvard University, 60 Garden Street, Cambridge, MA 02138; 

lars@cfa. harvard, edu 

ABSTRACT 

Recently, Pello et al. (2004) claimed to have discovered a galaxy at a redshift z = 10, 
and identified a feature in its spectrum with a hydrogen Lya emission line. If this 
identification is correct, we show that the intergalactic medium (IGM) around the galaxy 
must be significantly ionized; otherwise, the damping wing of Lya absorption by the 
neutral IGM would have greatly suppressed the emission line. We find either that the 
large-scale region surrounding this galaxy must have been largely reionized by z = 10 
(with a neutral fraction < 0.4) or that the stars within the galaxy must be massive 
(> 100 Mq), and hence capable of producing a sufficiently large H II region around it. 
We generalize these conclusions and derive the maximum Lya line flux for a given UV 
continuum flux of galaxies prior to the epoch of reionization. 
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1. Introduction 

Pello et al. (2004) have reported the discovery of a galaxy at a redshift z = 10, which is 
gravitationally magnified by a factor ~ 25 — 100 by the foreground galaxy cluster Abel 1835. The 
redshift identification is supported by the spectroscopic detection of an emission line at 1.338/mi, 
the photometric signature of an absorption trough at shorter wavelengths, and the lensing geometry. 
Here, we examine the far-reaching implications of associating this spectral feature with the Lya 
emission line from a galaxy at z = 10. In particular, we show that if hydrogen in the intergalactic 
medium (IGM) around this galaxy were neutral, then the damping wing of its Lya absorption 
would have greatly suppressed the Lya emission line of the galaxy. Our results provide a strong 
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incentive for obtaining higher quality data on this galaxy. More generally, we derive the maximum 
Lya line flux for a given UV continuum flux of galaxies prior to the epoch of reionization. Future 
detection of bright Lya lines from high-redshift galaxies can be used to improve the current lower 
limit on the redshift of reionization. 

Throughout, we adopt values for cosmological parameters derived from recent observations of 
the cosmic microwave background (Spergel et al. 2003; Seljak et al. 2003). For the contributions to 
the energy density, we assume ratios relative to the critical density of Q m = 0.27, = 0.73, and 
= 0.046, for matter, vacuum energy, and baryons, respectively. We also choose a Hubble constant 
Hq = 72 km s _1 Mpc _1 , and a primordial scale-invariant (n = 1) power-law power spectrum with 
ag = 0.9, where eg is the root-mean-square amplitude of mass fluctuations in spheres of radius 
8 h- 1 Mpc. 



2. Basic Absorption Parameters 

The absorption profile owing to H I in the IGM depends on the redshift z$ and halo mass M 
of the source, and parameters that characterize the production of ionizing photons by the galaxy. 
The first parameter, dN^/dt, is the total rate at which hydrogen ionizing photons from the galaxy 
enter the IGM. The second is the age of the source, ts, namely the period of time during which the 
source has been active (assuming a steady dN^/dt during that time, for simplicity). 

We consider ages in the range ts ~ 10 7 -10 8 yr, comparable to the dynamical time within the 
host galaxies of interest. For given values of M and ts, we find it useful to express dN^/dt as 
follows: 

,/Y M O, /V,... 

(1) 



dJV 7 _ M n b N ion 



dt m p Q m ts 

where m p is the proton mass, and iV"i on gives the overall number of ionizing photons per baryon 
in the galaxy (ignoring helium). To derive the expected range of possible values for N ion , we note 
that it is determined by 

N ion = N J Uf esc , (2) 

where we assume that baryons are incorporated into stars with an efficiency of /*, that iV 7 ionizing 
photons are produced per baryon in stars, and that the escape fraction (into the IGM) for the 
resulting ionizing radiation is / esc . When necessary, we adopt /* = 10% (consistent with a compar- 
ison of numerical simulations to the observed cosmic star formation rate at low redshift [Hernquist 
& Springel 2003]), and consider / esc in the full range 0-1. As explained below, our main conclusions 
are insensitive to our choices for the values of /* and / esc - 

The parameter JV 7 depends on the initial mass function (IMF) of the stars. In what follows, 
we consider two examples. The first is the locally-measured IMF of Scalo (1998), which yields 
[using Leitherer et al. (1999)] iV 7 = 4300 for a metallicity of 1/20 of the solar value. The second 
is an extreme Pop III IMF, assumed to consist entirely of zero metallicity M > 100 M & stars, 
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which yields iV 7 = 44000 based on the ionization rate per star, the stellar spectrum, and the main- 
sequence lifetime of these stars (Bromm, Kudritzki, &; Loeb 2001). Note that the star formation 
rate (hereafter SFR) corresponding to eq. (1) is 

— (Jk) (ife)" (&) «*- 1 



3. Ionization and the Damping Wing 

An ionizing source embedded within the neutral IGM can ionize a region of maximum physical 

size 

G^r m" mr (»r « 

assuming that recombinations are negligible and that all ionizing photons are absorbed by hydrogen 
atoms. In this Letter we conservatively adopt this maximum size, since a smaller region would only 
strengthen our main conclusions (for more elaborate discussions, see Haiman 2002; Santos 2003). 

When dealing with Lya emission and absorption, it is useful to associate a redshift z with a 
photon of observed wavelength A G bs according to A G bs = A a (l + z), where the Lya wavelength is 
X a = 1215. 67A. The optical depth rd am p owing to the IGM damping wing of neutral gas between 
Z\ and Z2 (where z\ < Z2 < z) can be calculated analytically (Miralda-Escude 1998). For a source 
above the reionization redshift, Z2 is the redshift corresponding to the blue edge of its H II region, 
and z\ is the reionization redshift. We define A A = X a (z — 22), which is related to the physical 
distance AD between the source and the H II region edge by 

aa \ A + z\~ 5/2 fn m h 2 ^- 1/2 



AD = - 1M {m){—) iwj Mpc - (5) 

In the limit where Z2 is very close to z while z\ is much farther away, the damping optical depth 
can be accurately approximated as 



1.16 Mpc fn b /n m \ f AA 11 

r damp - \^ w - ) |1- — - 1Q 3 A l - - log 

which is independent of z\. 



2.68 x 10 4 A 1 + 



AA 11 



(6) 



4. Lya-emitting Galaxies 

Next we discuss the use of Lya-emitting galaxies to probe reionization. To study the de- 
tectability of high-redshift galaxies, we must convert the SFR of a galaxy in a halo of a given 
mass [see eq. (3)] to a flux level for both the Lya line and the UV continuum. Assuming that 2/3 
of the ionizing photons that are reprocessed (i.e., that do not escape from the galaxy) lead to a 
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hydrogen atom in the 2p state and thus result in the emission of a Lya photon (Osterbrock 1974), 
the rest-frame conversion factor between the SFR and the luminosity of the Lya line is 



L(Lya) = 1.8 x 10 42 (1 - / esc ) ( ( SFR T ) erg s" 



(7) 



AT, 



where this conversion is valid for any ts > 10 7 yr (since the high-mass, short-lived stars dominate 
the ionizing intensity). Using the photon production rate in Leitherer et al. (1999) for the Scalo 
(1998) IMF, and the stellar spectrum and lifetime in Bromm et al. (2001) for the extreme Pop III 
IMF, we find that the continuum flux can be converted (in the rest frame) as follows, depending 
on ts and source rest-frame wavelength A: 

F U (X) = 10 27 7V cont ( SFR A erg s" 1 Hz" 1 , (8) 
V 1 M yr V 

where 

10. , for a Scalo (1998) IMF, t s = 10 s yr, A = 1500A , 

8.7 , for a Scalo (1998) IMF, t s = 10 8 yr, A = 2000A , 

5.5 , for a Scalo (1998) IMF, t s = 10 7 yr, A = 1500A , 

4.6 , for a Scalo (1998) IMF, t s = 10 7 yr, A = 2000A , ^ ' 
5.9 , for an extreme Pop III IMF, t s > 3 x 10 6 yr, A = 1500A , 

3.8 , for an extreme Pop III IMF, t s > 3 x 10 6 yr, A = 2000A . 

Note that for a given SFR and / esc , an extreme Pop III population produces a Lya line ~ 10 times 
more luminous than does a Scalo (1998) IMF, but the UV continuum is comparable for the two 
IMFs. The reason is that high-mass, short-lived stars strongly dominate the production of ionizing 
photons, but in the Scalo IMF case, lower-mass and longer-lived stars contribute substantially to 
the continuum. 

To evaluate the impact of absorption on the Lya line, we require an estimate of the width of 
the line. Assuming a Gaussian profile, we estimate the one-dimensional velocity dispersion a v by 
a v « V c / y/2, where V c is the virial circular velocity of the halo; this yields 

M y/3 /i + z y /2 fn mh 2\^ 



a v = 26.9 I — — km s" 1 . (10) 

\io 9 M y V 11 / V 0.14 1 



The observed scale-length <j\ is related to a v by 

ri) v 1 1 



rr, -22.1 ( ( i-— ) 1 kms 1 . (11) 



5. Results 



Identifying the redshifted Lya line based on its properties and a strong photometric continuum 
break, Pello et al. (2004) measured an amplified line flux for the candidate z = 10.0 galaxy of 
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(4.1 ± 0.5) x 1CT 18 erg cm -2 s _1 and a UV continuum flux density of (3.6+J' °) x 10~ 30 erg cm~ 2 
s _1 Hz -1 averaged over the 1360-1640A rest-frame band, and (2.3^{j y) x 10~ 30 erg cm -2 s _1 Hz -1 
averaged over the 1830-2090A rest-frame band. Correcting for the cluster lensing amplification 
factor of ~ 25, they estimated a SFR of 2-3 M & yr" 1 based on the continuum strength and 0.03- 
0.09 M yr~ x based on the line flux. This discrepancy, they suggested, indicates a factor of ~ 40 
absorption of the line, in principle owing to both H I absorption and dust extinction, but in this case 
they noted that the steep UV slope indicates negligible extinction. The observed line is somewhat 
resolved, allowing Pello et al. (2004) to place an upper limit of ~ 200 km/s on a v , with best fits 
obtained for a v < 60 km/s. 

We reconsider the interpretation of the observations of Pello et al. (2004), accounting for 
different possible IMFs and source ages. We find a luminosity distance to the galaxy of 105 Gpc, 
yielding a rest-frame luminosity (demagnified by a factor of 25) of 2.2 x 10 41 erg s" 1 in the line 
and a (demagnified) rest-frame luminosity density of 1.7 x 10 28 erg s -1 Hz -1 and 1.1 x 10 28 erg 
s _1 Hz -1 in the two UV continuum ranges, respectively (including the extra redshift factor for 
frequency bands). Estimating the overall conversion factor with the values for A = 1500A and for 
A = 2000A, respectively, we find using eq. (8) the SFR implied by the continuum flux, which then 
implies through eq. (7) an intrinsic Lya line luminosity (before absorption). We also estimate a 
halo mass using eq. (3), assuming = 0.1. 

Table 1 shows our derived parameters for the z = 10.0 galaxy and its host halo, depending 
on the IMF and the source lifetime. The Table gives the intrinsic line luminosity expected in 
the absence of H I absorption or dust extinction; also given are the expected intrinsic a v based on 
eq. (10), and the maximum physical size of the H II region based on eq. (4). Given the observations, 
the derived quantities given in the Table are independent of the star formation efficiency /*, except 
for the estimated halo mass and a v . Note that the observed ratio of the values of the continuum flux 
in the two measured wavelength regions is 1.6 ±0.7, compared to a predicted ratio of 1.2 for a Scalo 
IMF and 1.6 for a Pop III IMF. Thus, the steep UV slope is more consistent with a Pop III IMF, 
but the discrepancy with a Scalo IMF is not highly significant. In either case, this measurement 
suggests that dust extinction is not strong (compare discussion in Pello et al. 2004). For each case 
in the Table we average the SFR derived from the two UV continuum measurements. The errors in 
these two measurements introduce errors of ~ 20% in the predicted Lya line luminosity and ~ 7% 
in the predicted -R m ax- 

For a given value of / CSC) the rest-frame Lya line luminosity can be predicted from the UV 
continuum by accounting for the factor of 1 — / esc (which arises since ionizing photons produce 
Lya only if they are absorbed within the galaxy and reprocessed), and adding the damping wing 
absorption which depends on i? max (which depends on / csc as well). For each case in the Table, 
we vary / esc over all possible values and find the maximum possible Lya line luminosity L^. If 
we allow a 2-cr measurement error in the UV continuum flux, then the value of L^ ax becomes, 
respectively for the four cases shown in Table 1, 2.1 x 10 39 erg s _1 , 2.6 x 10 40 erg s~\ 1.4 x 10 42 
erg s _1 , and 9.3 x 10 42 erg s -1 . These numbers can be compared to the observed value of 2.2 x 10 41 
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erg s _1 (with a 2-a lower limit of 1.7 x 10 41 erg s -1 ). Thus, the Scalo (1998) IMF predicts a much 
weaker line than that observed and is ruled out, assuming that the IGM surrounding the H II region 
is neutral. An extreme Pop III IMF, however, is easily consistent with the observations even for a 
short burst (ts = 10 yr). 

Each line in the Table provides the mass of the source halo. Interestingly, numerical simulations 
indicate that the multiplicity function of star formation at z ~ 10 peaks at halo masses ~ 10 9 M Q 
(Springel k, Hernquist 2003), comparable to our inferred values. Using a halo mass function based 
on numerical simulations (Sheth & Tormen 1999, 2002), we find that the mean physical separation 
between halos at or above the inferred mass is 0.21 Mpc, 0.54 Mpc, 0.21 Mpc, and 0.82 Mpc, 
respectively, for the four cases in the Table. 

We have found that the Scalo (1998) IMF is inconsistent with the observations if the H II 
region is produced by the source galaxy and is surrounded by fully neutral IGM. There are two 
possible ways to resolve this inconsistency. One is to assume a larger surrounding H II region that 
is produced by the collective ionizing flux from nearby galaxies. The required minimum size for 
producing the observed line flux (with / esc = in order to maximize the intrinsic line flux) is 0.35 
Mpc (ts = 10 7 yr) or 0.43 Mpc (ts = 10 8 yr). Given the above halo abundances, such a sphere will 
typically contain a number of equally or more massive halos, i.e., an average number of 20 and 2, 
respectively. Barkana & Loeb (2004b) showed that the fluctuations in the number of halos in such a 
volume are much larger than Poisson fluctuations, owing to biased galaxy formation which greatly 
amplifies large-scale density perturbations. In particular, the typical 1-a fluctuation in the number 
of halos is a factor of 3 and 5 in these two cases, respectively; thus it is possible that the selection 
effect of the observability of the Lya line has revealed a galaxy in a region that had been previously 
reionized, e.g., by a positive 2-a fluctuation in the local number density of galaxies, while most of 
the rest of the universe remained neutral at that redshift. 

Another possibility is that the H II region has the size given in the Table, but the IGM beyond 
the H II region is on average only partially ionized and thus produces a weaker damping wing than 
we have assumed thus far. If as before we maximize the flux over all possible values of / C sc, then 



Table 1. Derived parameters for the Pello et al. (2004) galaxy (see text) 



IMF 


ts 

[yr] 


SFR 

[M e yr- 1 ] 


Halo mass 
[M & ] 


Lya luminosity 

[(1 - /esc) erg s" 1 ] 


[km/s] 


^max 

[/cVc 3 Mpc] 


t a 

max 

[erg s" 1 ] 


Scalo 


10 7 


2.8 


1.6 x 10 9 


5.0 x 10 42 


32 


0.16 


5.1 x 10 38 


Scalo 


10 8 


1.5 


8.8 x 10 9 


2.7 x 10 42 


56 


0.28 


9.4 x 10 39 


Pop III 


10 7 


2.9 


1.7 x 10 9 


5.3 x 10 43 


32 


0.36 


5.5 x 10 41 


Pop III 


10 8 


2.9 


1.7 x 10 10 


5.3 x 10 43 


69 


0.77 


4.8 x 10 42 



-7- 



-15 - 



-20 



-25 



-20 - 



i — i i 1 1 1 1 




I I I I I lll| 1 I I I I lll| 1 I I I I lll| 1 I I I 1 1 1 1 




JQ-32 10-31 1Q-30 

F„(1500 A) [erg/cm Vs/Hz] 





10-89 



Fig. 1. — Maximum observable flux of the Lya line versus UV continuum flux (at a rest-frame 
wavelength of 1500A) of a galaxy prior to reionization (i.e., assuming that the H II region of the 
galaxy is surrounded by a neutral IGM). We assume a narrow emission line (< 100 km s" 1 ). We 
consider a Scalo (1998) IMF (bottom panel) or a Pop III IMF (top panel), at redshifts 7 (short- 
dashed curves), 10 (solid curves), or 13 (long-dashed curves). In each case, a pair of curves are 
shown, corresponding to tg = 10 8 yr (higher) or t$ = 10 7 yr (lower). Also shown are the observed 
flux values for the Pello et al. (2004) z = 10 galaxy, demagnified by a factor of 25 (denoted by x). 
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the maximum allowed neutral fraction is 0.27 (t s = 10 7 yr) or 0.37 (t s = 10 8 yr). 

The damping wing of a neutral IGM limits the maximal flux of the Lya line from any galaxy 
prior to reionization. Figure 1 illustrates this upper limit as a function of the UV continuum flux of 
that galaxy, which affects the size of the H II region around it. We have maximized the line flux over 
/esc and calculated the damping wing at the center of the line assuming that the IGM is neutral and 
that its opacity does not change significantly across the width of the line. The latter assumption 
is well satisfied for lines with velocity widths < 100 km s _1 . Our flux limits are independent of the 
star formation efficiency /*. The galaxy discovered by Pello et al. (2004) (denoted by x) should be 
compared to the solid lines (indicating results for z = 10); its detected Lya line flux requires that 
the IGM be ionized for a Scalo (1998) IMF (bottom panel) but is consistent with a neutral IGM 
for a Pop III IMF (top panel). 

6. Discussion 

Figure 1 shows that if the galaxy discovered by Pello et al. (2004) indeed has a redshift z = 10, 
then either its stars are very massive (> 100M Q ) and hence capable of creating a large H II region 
around it, or the large-scale IGM around it has already been mostly reionized (with a neutral 
fraction < 0.4). The ratio of the galaxy's UV continuum flux in two bands is consistent with 
the first possibility, but better determination of the continuum spectrum or detection of other 
recombination lines (Bromm, Kudritzki, & Loeb 2001; Tumlinson, Giroux, & Shull 2001; Oh, 
Haiman, & Rees 2001) is required to make a convincing case. For example, a stellar population 
with a top-heavy IMF will produce substantial flux in He II Lya (303. 9A) and He II A1640. A 
search for features in the spectrum of a redshift z = 10 galaxy at A Q bs = 3343A and 1.8/xm would 
constrain the first interpretation. The discovery of other galaxies may test the second possibility, 
since the damping wing of a neutral IGM limits the maximum Lya line flux of any galaxy prior to 
reionization based on its UV continuum flux. 

The imprint of the damping wing is insensitive to a variety of complicating factors. Clumping 
of gas and its infall onto the galaxy have a minor effect, since the damping wing averages the IGM 
over large physical scales (~ 1 Mpc). Moreover, the characteristic size of the H II region around 
a bright galaxy at z ~ 10 is an order of magnitude larger than the scale of its infall region. The 
characteristic peculiar velocity of galaxies at z ~ 10 is an order of magnitude smaller than the 
Hubble velocity at the edge of their H II regions, and can only significantly change the effect of 
resonant absorption on the line profile. We have ignored resonant absorption from H I within the 
H II region since it would not suppress even the blue wing of the Lya emission line if the galaxy had 
a peculiar velocity (in the direction of recession) larger than the line width (thus allowing the line 
to appear nearly symmetric); such a peculiar velocity may not be improbable within small samples 
because it makes the line more easily observable and adds a selection bias. In any case, the addition 
of resonant absorption would only strengthen our limits on the Lya line flux before reionization. 
Dust extinction in the host galaxy would affect ionizing photons (and thus the Lya line flux) more 



-9- 



than it would the continuum at longer wavelengths. Based on the observed UV continuum, we 
would have predicted a weaker line with the inclusion of extinction, and so our argument would 
have been stronger. 

We have considered the lowest (most probable) value within the range of lensing magnification 
factors inferred by Pello et al. (2004). A higher value would again strengthen our limits. For 
example, a magnification factor of 100 would imply that the line was 4 times weaker intrinsically, 
but the predicted line from the UV continuum would then be weakened by more than a factor of 
4 since the smaller H II region would imply a damping wing that suppressed more of the line flux. 

Current observational constraints on reionization provide an inconsistent picture. On the one 
hand, the large-scale polarization anisotropics of the cosmic microwave background measured by 
WMAP imply a reionization redshift of 10-20 (Kogut et al. 2003), while on the other hand, the 
extent of the H II regions around the highest redshift quasars indicates 1 a significantly neutral IGM 
at z ~ 6.4 (Wyithe & Loeb 2004). Moreover, the IGM should have been cooler than observed at 
z ~ 3-4 if hydrogen had been fully reionized at z > 9 (Theuns et al. 2002; Hui h Haiman 2003). 
Theoretically, it is possible that the ionization fraction evolved in a complex, non-monotonic fashion 
owing to an early episode of Pop III star formation (Wyithe & Loeb 2003; Cen 2003). 

Existing observations do not enable us to discriminate between the two interpretations of a 
positive detection of a Lya emission line from the z = 10 galaxy of Pello et al. (2004). Sokasian et 
al. (2003a,b) have studied the process of reionization using cosmological simulations with radiative 
transfer; they have shown that for a Scalo-type IMF, reionization occurs over an extended redshift 
interval. For models in which the Lya optical depth at z ~ 6 matches the value inferred from the 
SDSS quasars, a substantial fraction of the mass in the universe (> 30%) is already completely 
ionized by z = 10. Since reionization proceeds "inside-out" affecting overdense regions first within 
the simulated volume (see also Gnedin 2000), we would expect the environment of the Pello et al. 
(2004) galaxy to be ionized. If the evolution had been more complex and included sources with 
population III stars prior to z = 12, an even larger fraction of the mass would have been ionized 
by z = 10 (Wyithe & Loeb 2003; Sokasian et al. 2003b). This picture could change, however, if 
reionization occurred "outside-in" (e.g., Miralda-Escude, Haehnelt <fe Rees 2000). As shown by 
Barkana & Loeb (2004), the fluctuations in the real universe should be far larger than indicated in 
existing simulations because of their limited box size. In reality, small fluctuations on scales larger 
than the box are greatly amplified through biased galaxy formation at high redshifts. 

The various scenarios for reionization can be tested with more data on Lya emitting galax- 
ies (supplementing recent work by Hu et al. 2002, 2004; Malhotra & Rhoads 2002; Kodaira 
et al. 2003; Santos et al. 2003; Kneib et al. 2004; Stanway et al. 2004; Barton et al. 2004), 
SDSS quasars (Fan et al. 2003), or gamma-ray burst afterglows (Loeb 2003; Barkana & Loeb 



1 The evidence for a possible Gunn-Peterson (1965) trough in these quasar spectra (White et al. 2003) is still being 
debated (Songaila 2004). 
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2004a). Over the next decade, it may also be possible to map directly the neutral hydrogen in the 
IGM through its 21cm line (see recent discussions by Zaldarriaga, Furlanetto, k Hernquist 2003; 
Morales k Hewitt 2003; Loeb k Zaldarriaga 2003; Gnedin k Shaver 2003; Sokasian, Furlanetto 
k Hernquist 2003) with forthcoming instruments such as LOFAR (http://www.lofar.org/) or SKA 
(http://www. skatelescope. org/) . 
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